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, nary systems up to the 5/2 post-Newtonian approximation of general rela- 

tivity (where radiation-reaction effects first appear). The approximate post- 
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Newtonian gravitational field might be used in the problem of initial condi- 
tions for the numerical evolution of binary black-hole space-times. On the 
other hand we recover the Damour-Deruelle 2.5PN equations of motion of 



compact binary systems. Our method is based on an expression of the post- 

H : 

Newtonian metric valid for general (continuous) fluids. We substitute into 

the fluid metric the standard stress-energy tensor appropriate for a system 

of two point-like particles. We remove systematically the infinite self- field of 

each particle by means of the Hadamard partie finie regularization. 
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I. INTRODUCTION 



The two purposes of the present paper are: 

(1) To obtain the gravitational field generated by a system of two point-like particles up 
to the so-called 5/2 post-Newtonian (2.5PN) order included, i.e. the order (t>/c) 5 where v 
denotes a typical value of the orbital velocity of the system. The 2.5PN field may be useful 
for setting up initial conditions for the numerical study of the coalescence of two compact 
(neutron stars or black holes) objects 

(2) To derive from the gravitational field the Damour-Deruelle |3|-[7| equations of motion 
of (compact) binary systems at the same 2.5PN order. Because the 2.5PN term in the 
equations of motion represents the dominant contribution of the radiation reaction force, the 
Damour-Deruelle equations play a crucial role in theoretically accounting for the decreasing 
of the orbital period of the binary pulsar PSR 1913+16 |9]^T2|. 

In addition, the present paper is motivated by the current development of the future 
gravitational-wave observatories LIGO and VIRGO. Specifically the aim is to derive with 
sufficient post-Newtonian precision the dynamics of inspiralling compact binaries (which 
are among the most interesting sources to be detected by LIGO and VIRGO). Numerous 
authors |T3"|-|18f have shown that the orbital phase of inspiralling compact binaries should 
be computed for applications in LIGO/VIRGO up to (at least) the 3PN relative order. 
Resolving this problem requires in particular the binary's equations of motion at 3PN order, 
since they permit to derive the binary's 3PN energy entering the left-hand-side of the energy 
balance equation on which rests the derivation of the phase. They are also needed for 
the computation of the 3PN gravitational flux entering the right-hand-side of the balance 
equation. Thus the 2.5PN equations of motion derived in 0-0] and in the present paper 
are not quite sufficient for the problem of inspiralling compact binaries, but the method we 
propose should permit to tackle in future work the problem of generalization to the next 
3PN order (see || for an attempt at solving this problem). 

The dynamics of a binary system of point-like particles modelling compact objects was 
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investigated successfully up to 2.5PN order by Damour, Deruelle and collaborators |3|-0, 
using basically a post-Minkowskian approximation scheme (i.e. G — > 0). 

In a first paper by Bel et al 0| (see also [[0J), the gravitational field and the equations of 
motion are obtained in algebraically closed-form to the second post-Minkowskian order (G 2 ): 
the field equations in harmonic coordinates are solved at first order by integrating the matter 
stress-energy tensor suitable to point-like sources (i.e. involving delta functions), and then 
the second-order gravitational field is constructed by iteration. The divergencies which arise 
due to the assumption of point-like particles are cured by means of a regularization process 
based on the Hadamard partie finie |20| (see |2l[] for an entry to mathematical literature). 
The equations of motion are obtained equivalently from the harmonicity condition to be 
satisfied by the metric, from the conservation of the (regularized) stress-energy tensor, or 
from the regularized geodesic equations. 

In sequential papers the post-Minkowskian equations of motion are developed up to 
the order G 2 /c 5 , i.e. neglecting any term of the order 1/c 6 when c —>■ oo and any term of the 
order G 3 when G — ■> 0. However, it is well known [ [2"2"| , |2"B| that in order to obtain the complete 



equations of motion to the dominant 2.5PN order of radiation reaction, the latter precision 
is not sufficient because of the occurence of terms coming from the third post-Minkowskian 
metric (G 3 ) which contribute to both 2PN (1/c 4 ) and 2.5PN (1/c 5 ) approximations. These 
terms of orders G 3 /c A and G 3 /c 5 have been added by Damour ||[7[], thereby completing the 
2.5PN binary's equations of motion. Let us refer to the above derivation of the dynamics of 
a binary system as the "post-Minkowskian" approach. 

When obtaining the cubic terms within the post-Minkowskian approach |],[F| the two 
objects are not described by standard delta functions but rather by a Riesz kernel [|4[] de- 



pending on some complex parameter A. For non-zero A this kernel has an infinite spatial 
extension, but reduces to the Dirac distribution when A —>■ 0. The metric is defined by com- 
plex analytic continuation from the A-dependent post-Minkowskian iteration. The physical 
equations of motion corresponding to point-like particles are obtained from the /1-dependent 
metric by taking the limit A — > at the end of the computation. It has been proved ||[7j 



that the limit exists up to the 2.5PN approximation (no poles ~ 1/A develop to this order 

The motion of two particles to the 2PN order (neglecting the radiation-reaction 2.5PN 
terms) is conservative, i.e. there exist ten integrals of motion corresponding to the Newtonian 
notions of energy, linear momentum, angular momentum and center of mass position. It has 
been shown [^,[7| that the constants of motion can be recovered by variation of a generalized 
Lagrangian depending on the positions, velocities and accelerations of the bodies (recall that 
generically, i.e. in most coordinate systems, a 2PN Lagrangian depends upon accelerations 
|27|). Adding up the radiation- reaction terms, one finds that the previously obtained binary's 
2PN energy decreases with time, and that there is quantitative agreement with the standard 
quadrupole formula |]2"8|]7|1 and with the observations of the binary pulsar ||9|-[T2|. 

Moreover, there have been two other lines of work which led to the complete 2.5PN 
dynamics of binary systems. One of these alternative approaches is based on the canonical 
Hamiltonian formulation of general relativity and the model of point-like sources. Such 
"Hamiltonian" approach was developed at the 2PN level in early works P§[ , PU| (see also 



but completely understood only later ]27[ (see [[33] for a review). Schafer 
completed the Hamiltonian approach to include the 2.5PN radiation-reaction terms (more 
recently the 3.5PN radiation-reaction terms have also been worked out f36j). The other 



method gives up the model of point-like sources and assumes from the start that the two 
bodies are extended, spherically symmetric and made of perfect fluid. Within such an 
"extended-body" approach the 2.5PN equations of motion were found |37|j38|] to be the same 
as obtained within the other treatments dealing with point-like particles (in particular the 
equations depend only on the two masses of the bodies, but not on their internal structure 
nor compactness). 

In the present paper, we add what constitutes essentially a fourth approach to the deriva- 
tion of the 2.5PN motion of binary systems, which can be qualified as "post-Newtonian", 
in contrast with the post-Minkowskian, Hamiltonian and extended-body approaches. With 
respect to the post-Minkowskian approach [§-0] we have essentially two differences: 



(1) Instead of implementing a post-Minkowskian algorithm to the third order and 
performing afterwards the post-Newtonian re-expansion, we start directly from a post- 
Newtonian metric developed to 2.5PN order and which is valid for any continuous matter 
stress-energy distribution ("fluid"). Note however that our post-Newtonian metric is defined 
in terms of retarded (Minkowskian) potentials, and most importantly matches a far-zone 
metric satisfying the correct boundary conditions at infinity, in particular the no-incoming 
radiation condition . 



(2) Instead of assuming a fictituous stress-energy tensor defined by means of analytic 
continuation using the Riesz kernel and letting the analytic- continuation factor going to 
zero at the end of the computation ||[7|], we substitute directly into the "fluid" metric 
the stress-energy tensor of point-like particles. This entails divergencies which are cured 
systematically by means of the Hadamard regularization pO| , |2T| (in this respect we follow 



By implementing the post-Newtonian approach, we recover the Damour-Deruelle equa- 
tions of motion Q-|7[|. To the investigated order we find that the post-Newtonian approach 
is well-defined and rather systematic. 

About the gravitational field (metric) generated by the two particles, we obtain an alge- 



braically closed-form valid everywhere in space-time up to the 2.5PN order p0| . Indeed, the 



most difficult term present in the metric at this order is a cubically non-linear term which 



can be explicitly evaluated pH] , f42|| . This yields the other motivation of the present paper, 
namely to provide the metric coefficients at 2.5PN order (in harmonic coordinates) in the 
form of some explicit, fully reduced functionals of the positions and (coordinate) velocities 
of the two masses. Let us point out that, very likely, the possibility of writing such a closed- 
form expression of the metric breaks down at the next 3PN order, where there remain some 
Poisson-type integrals which probably cannot be expressed in terms of simple functions. 

The plan of the paper is as follows. We start in Section II with the expression (derived in 
Appendix A) of the 2.5PN metric valid for general fluid systems. In Section III we explain 
our method for applying the fluid metric to the case of point-like particles. The metric 



potentials involve three types of terms which are evaluated respectively in Sections IV, V 
and VI (the most difficult, cubic, term being obtained in Section VI). The results for the 
potentials are relegated to Appendix B. In Section VII we present our expression for the 
binary's gravitational field. In Section VIII we finally obtain the (Damour-Deruelle) binary's 
equations of motion. 



II. THE 2.5PN METRIC FOR GENERAL FLUID SYSTEMS 

At the basis of our investigation is the expression of the metric generated by an arbitrary 
matter distribution described by the stress-energy tensor T^ v [f|3| . We assume that T^ u has a 
spatially compact support, and physically corresponds to a slowly- moving, weakly-stressed 
and self-gravitating system, in the sense respectively that |X°yx 00 | ~ g, \T l i /T m \ ~ e 2 
and U/c 2 ~ e 2 , where U denotes the Newtonian potential and e represents a small post- 
Newtonian parameter going to zero when the speed of light tends to infinity (e ~ 1/c). 
Throughout this paper we denote a post-Newtonian term of order e n by means of the short- 
hand O(n). 



Following |44j] it is convenient to define a mass density a which agrees in the case of 
stationary systems with the Tolman mass density to 1PN order. As it turns out, introducing 
such a mass density (and in addition the associated retarded potential) permits to formulate 
the 2.5PN metric in a rather simple fashion. Defining also some current and stress densities 
we pose 

^00 I rj^ii 

o = - 2 , (2.1a) 

c z 

(Ti= — , (2.1b) 
C 

o u - •/"' . (2.1c) 

The covariant conservation of the matter stress-energy tensor (V u T fJ-u = 0) entails the equa- 
tions of motion and continuity, which read (with relative 1PN precision in the equation of 
continuity but only Newtonian precision in the equation of motion) 



d t a + dm = \ (d t a u - ad t U) + 0(4) , (2.2a) 



c 2 



d t a t + d^y = adiU + 0(2) , (2.2b) 

where U is given by the standard Poisson integral: U = A _1 {— AirGa}. 

Actually, it is advantageous to use rather than the instantaneous potential U, a corre- 
sponding retarded potential V given by the retarded integral of the same source a: 

V{x, t) = n R l {-AnGa} = G [ — <r(z, t - |x - z|/c) . (2.3) 

To Newtonian order we have V = U + 0(2) ||45| . Similarly let us introduce the following 
other retarded potentials [fHf] : 

Vi = a R l {-AnGa.i} , (2.4a) 

Wij = n R l {-47rO((7y - 6 tJ a kk ) - d.Vd.V} , (2.4b) 

Ri = ^ {-47rG(ya l - V { a) - 2d k VdiV k - ^Vw} , (2.4c) 
X = n^^-ArrGVau + 2Vid t diV + Vd 2 t V 

+^(d t V) 2 - 28^0^ + W^v} . (2.4d) 

In addition, we shall often consider the trace of the potential W^, i.e. 

W ti = a,, 1 {torGau - d.Vd.V} . (2.4e) 

We are now able to express the usual covariant metric g^ v in terms of these retarded 
potentials to order 2.5PN, by which we mean neglecting all the terms of order 0(8) in g 00 , 
0(7) in goi and 0(6) in g^. We impose the harmonic or De Donder coordinate conditions, 
i.e. d v [y/—gg ,J ' u ] = 0, where g and g^ v are the determinant and the inverse of the matrix 
g^y. Actually, since we are working with an approximate post-Newtonian metric, the har- 
monic conditions need only to be satisfied approximately. To the 2.5PN order we have 
dA^9° u ] = 0(7) and d v W=gg*\ = 0(6). 

With these definitions the 2.5PN metric in harmonic coordinates takes the form: 
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9oo 



2 2 „ 
-1 + —V-—V 2 + 



c 2 c 4 <■" 



6 



4 



!h, = - — \ i --^R i + 0{f) , 
g l3 = ( 1 + hi + - 4 V 2 ) + -W l3 + 0(6) . 



X i;V, 4 — -!- 0(S) . (2.5a) 

(2.5b) 
(2.5c) 

For the sake of completeness, this rather simple result is proved in Appendix A. The sim- 
plicity in the formulation is due to our introduction of the mass density a as well as the use 
of retarded potentials H4^j39f1 . 

In the form ( |2.5| ) the metric contains only "even" terms explicitly (using the standard 
post-Newtonian terminology), that are terms with even powers of 1/c in goo an d <7y, and 
odd powers in g^. Indeed the "odd" terms, which are responsible for radiation reaction 
forces, are all hidden into the definitions of the retarded potentials (|2.3|)-(p^), and can be 
made explicit by expanding the retarded arguments with Taylor's formula. It is important 
in this respect to recall from [39] that ( [2.5D comes from the post-Newtonian expansion (valid 
only in the near zone) of some radiative metric defined globally in space-time and satisfying 
the no-incoming radiation condition at past null infinity. Hence the "odd" terms in the 
post-Newtonian metric (|2.5| ) correspond physically to the radiation reaction forces acting on 
an isolated system (with no source located at infinity). 

At 2.5PN order the harmonic-coordinate conditions are equivalent to the following dif- 
ferential identities: 



OA) ■■<)■, {Wij-^Wkk 



dAV i +- 



Ri + Wi 



0(4) 



0(2). 



(2.6a) 
(2.6b) 



These relations are in turn equivalent to the 1PN continuity equation and Newtonian equa- 
tion of motion given by ( |2.2j ). 

The potentials V and Vi are generated by the compact-supported source densities a and 
<7j. Similarly W%j and Ri involve a part generated by a compact-supported source, but also a 
part whose source is a sum of quadratic products of potentials V or Vi and their space-time 
derivatives. We shall refer to the former part of Wij and Ri as the compact ("C") part, and 



S 



to the latter as the "dVdV" or, sometimes, "quadratic" part. As for X, it consists of C and 
dVdV parts like for W^ and R iy but it also contains a term of different structure, generated 
by the product of Wij and dfjV [last term in fl2.4d|) l. This term itself can be split into two 



contributions arising respectively from the C and dVdV parts of W^. Since the C part of 
is a compact-supported potential similar to V and Vj, the corresponding term in X has 
actually the same structure as a dVdV potential. On the other hand, the dVdV part of Wij 
generates an intrinsically more complicated term in X we shall refer to as the non-compact 
("NC") or "cubic" part of the potential X. Precisely, our definitions are 

Wtj = wtf> + Mf V9V) , (2.7a) 
Rt = RP + R? VdV) , (2.7b) 
X = X™ + X^ 9 ^ + X( N °) . (2.7c) 

The compact parts are linear and quadratic functionals of the matter variables fl2.1|) . They 
read as 

W<>f ] = n- 1 {-4ttG(^ - 5 l3 a kk )} , (2.8a) 
RP = D^ 1 {-47rG(V<7< - V l( y)} , (2.8b) 
XW = n-^ {-4irGV* u } . (2.8c) 

The dVdV or "quadratic" parts involve both quadratic and cubic contributions. They are 

Wg V9V) = n R l {-dtVdjV} , (2.9a) 

Rf yay ) = n-j {-to - \d t vd-y} , (2.9b) 

Xovav) = n -i | 2y . dAV + y d 2y + 3 (fty)a 

-2d i Vjd j V i + w£ ) dlv} . (2.9c) 

Finally the only Non-Compact part is a cubic functional given by 

X^ = U R ^{wif VdV) dlv} . (2.10) 
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In practice the latter term is the most delicate to evaluate. Our terminology is slightly 
improper, as the so-called dVdV or quadratic potentials are as well as the NC potential 
generated by non-compact supported sources, and involve some contributions which are 
actually cubic in the matter variables. 



III. APPLICATION TO POINT-LIKE PARTICLES 



To apply the general 2.5PN metric presented in the previous section to the case of a 
point-mass binary we use the matter stress-energy tensor: 



T^(x, t) = vn(tX(tX{t)5(x - yi(t)) + 1 



(3.1) 



In our notation the symbol 1 <-> 2 means the same term but with the labels 1 and 2 ex- 
changed; 5 denotes the three-dimensional Dirac distribution; the trajectories of the two 
masses (in harmonic coordinates) are denoted by yi(t) and y2(t); the two coordinate veloc- 
ities are Vj(t) = dyi(t)/dt, v 2 (£) = dy 2 (t)/dt, and v% = (c, Vi), v% = (c, v 2 ); ii\ represents 
an effective time-dependent mass of body 1 defined by 

\ 



mi 



\ 



'99 pa' 



P a 
ViVl 



(3.2) 



c 2 ) x 

mi being the (constant) Schwarzschild mass, with g pa the metric and g its determinant E7 
Another useful notation is 

,2 



#l(t) = Pitt) 



1 + ^ 



(3.3) 



where v\ = v 2 . Both \i\ and fl% reduce to the Schwarszchild mass at Newtonian order: 
fix = mi + 0(2) and p,% = ni\ + 0(2). Then the mass, current and stress densities Q2-lp for 
two particles read 



a = p,i8(x - yi) + 1 <-> 2 , 
(Ti = Hiv\5(x - yi) + 1 2 , 
Oij = ^iv\v{5(y. - y x ) + 1 2 



(3.4a) 
(3.4b) 
(3.4c) 
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The stress-energy tensor of point-masses depends on the values of the metric coefficients 
at the very location of the particles. However the metric coefficients there become infinite 
and, consequently, we must supplement the model of stress-energy tensor (|3.1| ) by a prescrip- 
tion for giving a sense to the notion of the field sitting on the particle. In other words, we 
need a regularization procedure in order to remove the infinite self-field of point-like sources. 
The choice of one or another regularization procedure represents (a priori) an integral part 
of the choice of physical model for describing the particles. In the present paper we shall 
employ the Hadamard regularization based on the partie finie of functions admitting 

a special ("tempered") type of singularity. For discussion and justification of the use of 
the Hadamard regularization in the context of equations of motion in general relativity see 



Let us consider the class of functions F depending on the field point x as well as on two 
source points yi and y2, and admitting when the field point approaches one of the source 
points (ri = |x — yi| — > for instance) an expansion of the type 

F(x;y 1} y 2 ) = £ r^(n 1 ;y 1 ,y 2 )+0(r 1 ) (3.5) 

-k <k<0 

(where k e Z). We define the value of the function F at the source point 1 (and similarly 
at the source point 2) to be the so-called Hadamard partie finie, which is the average, with 
respect to the direction ni = (x — y-\)/r\ of approach to point 1, of the term with zeroth 
power of r% in (3.5). Namely [fi9[| , 

(F), = F(y i; yi , y 2 ) = J ^j^/o(n x ; yi, y 2 ) . (3.6) 

Furthermore, we use the Hadamard partie finie to give a sense to the spatial integral of the 
product of F and the Dirac delta function at point 1 (since F is singular on the support of 
the Dirac function). Indeed, we define 

J d 3 xF(x; y a , y 2 )5(x - Yl ) = (F)x , (3.7) 

where (F)i is given by (|3.6|) . 

11 



As a (trivial) example of the use of the Hadamard regularization, consider the potentials 
V and Vi to Newtonian order, given by 

y=^ + 0(2) + l<->2, (3.8a) 
^ = ^p-vl + 0(2) + (3.8b) 

They are infinite at point 1, but after applying the rule ( |3.6| ) we find 

Q0i = — + O{2) , (3.9a) 
(Vt)i = — 4 + 0(2) , (3.9b) 

where r 12 = |yi — y 2 | is the distance between the particles pOfl . Of course (V)i agrees with 
the standard Newtonian result. Applying the rule ( |3.7| ) we have, for instance, 

1 f 3 TA Gmim2 



- d 6 *aV = 1 z + 0(2) , (3.10) 
2 7 

also in agreement with the Newtonian result. 

We shall derive the binary's equations of motion in the so-called order-reduced form, by 
which we mean that in the final result all accelerations (and time-derivatives of accelerations) 
are replaced consistently with the approximation by the explicit functionals of the positions 
and velocities as given by the (lower-order) equations. So, in order to derive the 2.5PN 
equations of motion (and also the metric), we use the less accurate 1.5PN equations, given 
in harmonic coordinates by 

dv{ Gm 2 



dt ~ ni2 \ X "c 2 



Gmi .Gmo 9 9 , s 3. 

-5 - 4 2 - + v\ + 2v 2 2 - 4 ( Vl v 2 ) - -{n l2 v 2 f 

ri2 r i2 - 2 



12 

+ ^¥^12 [4(^0 - 3(n 12 ^ 2 )] + 0(4) , (3.11) 

c r 12 

with ni 2 = |yi — y 2 |/ri 2 and v i2 = vi — v 2 ; scalar products are denoted with parenthesis, 
e.g. (ni 2 t>i) = ni 2 .vi. The acceleration of body 2 is obtained by exchanging the labels 1 
and 2 (remembering that n 12 and v 12 change sign in this operation). 
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IV. THE COMPACT PARTS OF POTENTIALS 



In this section we derive the compact-supported potentials V and V{, and the compact- 
supported parts of the other potentials, W-^\ R\ C ^ and X^ c ^ defined by Q2.8|), for a binary 
system described by the stress-energy tensor (|3.1|) and the regularization (3.6). We need 
V to relative 2.5PN order, Vi to 1.5PN order, and the other compact potentials to 0.5PN 
order only. The main task is the computation of V, to which we focus mainly our attention. 
By Taylor-expanding at 2.5PN order the retardation inside the integral (|2.3| ) and using the 
mass density a in the form ( |3.4a|) , we get 

V = G (— - -d t (h) + A^(Am) - 



+ 



1 



2c 2 



6 c 3 



^dt(hrt)} + 0(6) + 1-2. 



(4.1) 



24c 4 tv ^'^ 120c 5 ' 
We recall that the effective mass jli given by ( |3.2|) - (|3.3|) is a function of time only. 

We start by deriving fii to 2.5PN order. Inserting the metric coefficients (|2.5|) into the 
expressions ( |3. 2|) - ( p73|) , we obtain 



fix = m x \\ + — 



+ 



-2(W«)i + \{V\ + \iVyvl - A{V t ) lV \ + 7 -vt 



+ 0(6) 



(4.2) 



where all the potentials are to be evaluated at the location of body 1, using the rule ( |3.6| ). 
We proceed iteratively. The first step consists in inserting into ( ^4.2| ) the potential V at body 
1 to Newtonian order (or, rather, 0.5PN order), which is simply the Newtonian result ( p.9a|) . 
This yields fii to 1.5PN order: 



fit = m 1 <{ 1 + — 



Gm,2 



+ 0(4) . 



(4.3) 



Since the time-derivative of /ii starts at 1PN order, namely 

Gm 1 7TT,2 



c 2 rf 2 



-2(n 12 ui) - (n 12 v 2 )\ + 0(4) , 



(4.4) 



we see that the first odd power of 1/c in V arises at 1.5PN order. Furthermore, using the 
constancy of the center of mass velocity, one can check that the first odd term in the gradient 
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of V arises at 2.5PN order (it contributes to the dominant radiation reaction force). From 
fl4.ip and ( |4.3| ) we deduce the value of (V)i up to 1.5PN order: 

1 , 



Gm 2 f 1 
" + - 



+ 



ri2 L c 
4G 2 mim2 



3Gmi 2 

.—-+2v 2 -- { n l2 v 2) 



3c 3 r 2 



(7112V12) + 0(4) . 



(4.5) 
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In addition to (V)i, we need (V^)i already given by (|3.9b|) , and the value at point 1 of the 
trace of Wij to 0.5PN order. The trace Wn is much simpler than the potential itself, and 
from (|2.4ej ) we derive the expression 



(['„ = A -1 \ 8nG [an - -aV 



2 



Under this form all integrals are compact-supported; at this order, we can insert V = 



U + 0(2). The odd term 0(1) is a mere function of time. From ( 4~q) we get immediately 



nil 



Grri2 
r\2 



Gm\ Gm 2 



- 2vt 



1G 2 m\m 2 



cr\ 2 



\n\2Vvi) 



0(2). 



T12 2r i2 

The effective mass [L\ at 2.5PN order is readily obtained from the previous relations: 

Gm 2 3 



(4.7) 



Aii = mi jl + J 



+- 



2 

Gm 

Tl2 



T\2 

1 

2 Vl 



2 1 



> / 1 .-, , 1 . , 9 1 Grrii 3 Gm 2 

-' " 4{v 1 v 2 ) + 2v 2 2 + -{n l2 v 2 f - 



+ 



2 ri 2 2 na 



8G 2 mim 2 



3c 5 r 2 
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niavia) + 0(6) 



(4.8) 



from which we straightforwardly deduce V to 2.5PN order. The only point is to compute 
the numerous (up to five) time-derivatives of r\. This gives rise to many terms depending 
on the acceleration and its time- derivatives, that we reduce order by order by means of the 
binary's 1.5PN equations of motion (since an acceleration already arises in the 1PN term of 
V). Once fully reduced, the result for V is still lengthy so we relegate it (together with all 
the relevant results for the potentials) in Appendix B. 

The potential Vi to 1.5PN order, and the other compact potentials Wjj , and 
to 0.5PN order, are obtained in the same way. As an example we give 
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w, 



(C) 



Grrii 
ri 



(y\v{ - 5^) 



+ 



171x1712 



crf 2 



(* j) 

"12^12 



- ^'(niavaa) + 0(2) + 1 <-> 2 . (4.9) 



V. THE QUADRATIC PARTS OF POTENTIALS 

By the definition ( [2.9Q the quadratic or dVdV potentials have their sources made of 
quadratic products of (derivatives of) the compact-supported potentials V, Vi and W}j. 
All the dVdV potentials are to be computed to 0.5PN order, which means in particular 
that we can replace in the sources V and V{ by the Newtonian-like potentials U and Ui 
(but we must beware of the fact that Wffl given by fl4.9|) involves a 1/c correction). For 
all the dVdV potentials we proceed similarly. We work out the sources using ( p.8|) and 
(4.9) and obtain some "self" terms, proportional to m\ and m 2 , together with some "inter- 



action" terms, proportional to m\m 2 . Time-derivatives are changed to spatial derivatives 
thanks to d t (l/ri) = v\du(l/ri) and df(l/ri) = a\du(l/ri) + v\v{duj (1/ri), a\ denoting 
the acceleration and du the partial derivative with respect to y\. In the interaction terms, 
we leave the spatial derivatives un-expanded, whereas, in the self terms, they are developed 
and "factorized" out in front of the terms. In the latter operation, we should remember 
that within the standard distribution theory the second spatial derivative of \ jr\ involves a 
distributional term |5I|,ETIHSI1 : 



A(i) =-47T(J(x- yi ) , (5.1a) 
d i{-) =- L j* y^(x- yi ) . (5.1b) 



Two examples of such a treatment of sources are 
G 2 m\ / 2 . ■ \ / 1 \ 



+ G ,2 m 1 m 2 9 H 9 2j ( ) + 0(2) + 1^2, (5.2a) 



r±r 2 



VdfV = G l m\ 



/ 1 



+ G 2 m 1 m 2 (aldu + v^d 2 Uj ) [—- J + 0(2) + 1^2. (5.2b) 
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We apply the Poisson integral on the source terms treated in the previous manner. 
Consider firstly the distributional terms, such as the one in the self part of V&*V. Although 
this term is ill-defined, because involving the product of the Dirac distribution 5(x — yi) by 
1/ri which is singular when x — * yi, the Poisson integral is computed unambiguously with 
the help of the Hadamard regularization ( |3.7| ), yielding zero in this case: 

dh 

z||z - yi| 



A 



-i 



47T 



J |x — z| |z — Vi | 

= . (5.3) 



1 



x - z z - yi 



z=yi 

For the computation of all non-distributional terms in the dVdV potentials, we take the 
example of 



W i3V9V) = n _! { _ 9iVdjV} 

= A- 1 {-diVdjV} + -^-^ J d 3 x {-diVdjV} + 0(2) , (5.4) 



whose "source" is given by ( 5.2a ). The Poisson integral of the self-terms can be readily 
deduced from A(lnrx) = 1/rf; on the other hand, that of the interaction terms is obtained 
by solving the elementary Poisson equation 

Ag = — . (5.5) 
r x r 2 



The solution is known [52] : 



g = \nS; S = r 1 + r 2 + r 12 . (5.6) 



The computation of the 1 / c-term in ( |5.4| ) involves essentially the spatial integral of l/r^ . 
Since it is divergent due to the bound at infinity (i.e. when r = |x| — > oo), we first compute 
the finite integral defined by integration over a ball of constant finite radius TZ. By writing 
the integrand as l/r^ = Ag and using the Gauss theorem, we transform the integral into 
a surface integral over the sphere of radius 71, 



f f± = j d 3 x Ag = I dQ (r 2 d r g) , (5.7) 

J|x|<7^ Tir 2 J\yL\<n Jr=1l 
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with d r = n % di. Into the latter surface integral we can replace the function g by its expansion 
at infinity computed from ( |5.6| ): g = ln(2r) + [— (ny\) — (ny 2 ) +r 12 ]/2r + 0(l/r 2 ). Neglecting 
the terms which die out in the limit TZ —>■ oo, we get 

1 ' ** . K + Iil + o(l-) . (5.8) 



47T J\x\<TZ T\T2 2 \TZ, 

As we can see, the divergent part of the integral is simply a constant, which will therefore 
vanish after application of the spatial derivatives dud 2 j in front of the term. This shows 
that we are allowed to use in this computation the finite coefficient in the right-hand-side of 
Q5.8D , which is nothing but the partie finie in the sense of Hadamard of the initially divergent 
integral f53|| . Setting yi = y 2 in (|5.8j ), we infer that the divergent integral of 1/rf can be 



replaced by zero. Then, (|5.8|) together with the last fact lead to 

_ _L J rf 3 x d . vd .y = G 2 m im2 ^ _ ^ + 0(2) (5 g) 

Gathering those results, we thereby obtain the looked-for potential as 



rtr(avdv) G 2 m\ f 2 5 ij \ 2 

W>j = — y In n + — -G m x m 2 ^ 



2^, 

(i j) J- / ij 

n\ 2 v 12 - 



G 2 m\m 2 



cr\ 2 



\ (3n? 2 - (ni2«i 2 ) 



+ 0(2) + 1^2, (5.10) 



where we have ^ = dud 2 jg. The first two terms are in agreement with a result of [J2 
All the (9y(9^-potentials are calculated in this way to obtain the complete expressions of 
potentials presented in Appendix B. 

Ending this section we list some formulas which are useful in the derivation of the dVdV- 
potentials, and even elsewhere. The first-order spatial derivatives of g read 

i9 = d lt g = - w l + n " , (5.11a) 
9i = d 2l g = -"S- n " , (5.11b) 
Oig = -ig -9i = — g — , (5.11c) 

and second-order spatial derivatives are (with n^ 2 = n\ 2 n\ 2 ) 
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_ , n\ J 2 - 5 ij n l l - 5 ij (n\ 2 - n[)(n{ 2 - n{) , uo , 

,9 = cr Uj g - — — - , (5.12a) 



= ^ = —^5 & ' (5 ' 12b) 

a = a a , ' (S ' J i ( ra i2- 7l i)( ra i2 + "2) r , 12 x 

iflfj = Ch^P - 1 — 2 • [b.iZC) 

Contracting with the Kronecker 5^, the relations (|5.12 ) become ug = A\g = l/riri2, gu = 
A 2 g = l/r 2 ri2, and, more interestingly, 

i9i = \(- " —) ■ (5-13) 



5 lj n 2 — 5 tj (n\ 2 + n 2 )(n{ 2 + n 2 ) 



1 + {n 1 n 2 ) 


7"l + T 2 - 


ri2 


S 


2r x r 2 




1 - (nin 12 ) 


_ n + r 12 - 


■ r 2 


5 


2riri2 




1 + (n 2 n 12 ) 


r 2 + r 12 - 


■ n 


5 


2r 2 ri2 





2 \r 1 r 2 r 2 r 12 r x r 12 
This simple result is a straightforward consequence of the helpful formulas 

(5.14a) 
(5.14b) 
(5.14c) 

Finally, we find the values of g and its derivatives at the location of body 1 (say) according 
to the Hadamard regularization (3.6): (g)\ = ln(2ri 2 ), together with 

Ugh = tt^ ; (9ih = -— , (5.15a) 

-5^ + 2ng 1 
Mjh = w~2 5 WWi = -7J-2- > (5.15b) 

z 'l2 z 'l2 

(tfS)i = —^J2 5 \9ijh = - 2 • (5.15c) 

^'12 '12 

These formulas are extensively used when getting the potentials at 1 (see Appendix B). 



VI. THE NON-COMPACT POTENTIAL 



The so-called non-compact potential is defined by ( |2.10| ) as the retarded integral of a 
source composed of the product of a V-type potential and W^ V9V \ the latter potential 
being itself given as the retarded integral of a source made of a quadratic product of Vs. 
Due to this purely cubic structure, one would expect a priori that the computation of the 
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non-compact term represents a non-trivial task, and even that it is not at all guaranteed 
that this term can be expressible with the help of simple algebraic functions (algebraically 
closed- form) . Rather surprisingly, the NC potential turns out to accept an algebraically 
closed-form up to 0.5PN order. As a result, one can find its explicit expression, valid for 
any source point x (the value when the source point sits on a particle yi^ following from 
the regularization process). To Newtonian order the closed-form expression of the NC term 
has already been obtained in [f41~l , [42| by combining some technical results derived earlier in 
|54] , |3llj55| . We shall present here a slightly different but totally equivalent form of X( NC ) 
at Newtonian order, and add to this the 0.5PN correction. Very likely the 1PN and higher 
corrections in X( NC ) do not admit any algebraically closed-form all over space-time, but the 
regularized values at the location of the two bodies can probably be carried out explicitly 
(these values are needed when investigating the equations of motion). 
To 0.5PN order the non-compact potential reads as 

£ (NC) = n -l {wff V " V %V} 

= A- {wr SV %Y} + l± I rfSc {W^ V %V} + 0(2) , (6.1) 

T (dVdV) 



where to this order V can be replaced by ( |3.8a| ) and by ( |5.10| ). The cubic source is 



easily obtained thanks to ( |3.8a| ) and ( |5.10| ). Hence we arrive at 



w [9Vdv) d 2 v = 



3 + ~2 6 (*-y^ 



G 3 m?m2 ll s ±^A _ lag (I) al tori _ 2 9 i (I) , 9 ; 



G 3 mjm 2 [ i ■ 1, ij - a /2 

12 L ^ J V ' 1 , 

+ 0(2) + 1^2. (6.2) 



We compute the Poisson integral of ( |6.2| ). The (ill-defined) distributional term in the self 
part (ex m\) of (|6.2j ) is treated unambiguously using the rule ( |3.7| ) and does not contribute 
to the Poisson integral. On the contrary the distributional term in the interaction part 
(ex m\m2) is well-defined and gives a net contribution. Easy terms are obtained from the 
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facts that A(l/rf) = 6/rf and A(ri) = 2/ri. The difficult point is to find the solutions of 
the two Poisson equations 



^1 = 2^.(1)^ 
AH, = 2 d% (i) ^ . 



(6.3a) 
(6.3b) 



We use the same notation as in |g2[ except that we add a subscript 1 to distinguish a function 
from its image obtained by the exchange of bodies 1 and 2. Remarkably, the solutions of 
( |6.3|) can be written down everywhere in space-time under explicit form f41||4"2"| ,56 : 



-A- Ai 
2 



In r! 



r 2 



+ -A 2 



#i = -Ai 
2 



lnr 12 

?~2 



+ 



r 2 



2r\ 2 r\ 



+ 



rf 2 r 2 



(6.4a) 



+ did 2i 



q lnri 
— + - Ai 

n ria 



lnri2 lnri 



-0 



2r 



12 



— d 2i [{dig)x] - 



r-2 



T\ " ' " ' ' 2r\r\ 2 

They are equivalent to the expressions given by the equations (3.48)-(3.49) in [ 4^| . By 
expanding all derivatives, we come to the completely developed forms 



(6.4b) 



1 



1 



K x = — 3 + , 
r 2 r 2 rf 2 



r\r 2 



+ 



2r?r 2 



+ 



12 



1' 12 

r 2 



2 r 2 r 3 



+ 



+ 



2 

r 2 



2r? 



2 ' 



2 r 12 



+ 



3rf 



+ 



2rf 4r{ 2 4r'(ri 2 2r'{r'{ 2 ' 2rfr i2 ' 4r^rf 2 ' 2r^r^ 2 2r\r\ 2 
(K 2 and i7 2 are obtained by exchanging ri and r 2 in the right-hand-sides). With the solutions 
( |6.4| )-( |675D , we control the NC potential at the Newtonian approximation. 

Next, we compute the spatial integral of (|6.2| ) entering the 0.5PN correction in the NC 
potential. We must evaluate essentially the spatial integrals of the two source terms in the 
right-hand-sides of (|6.3|) . We proceed like for the integral of \jr\r 2 in (|5.7|) - (|5.8|) . Namely, 
we integrate over a ball of constant radius TZ, and use the function K\ to transform the 
volume integral into a surface integral over the sphere r = TZ: 



(6.5a) 
(6.5b) 



(6.6) 
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From the developed expression of K\ given by ( p.5a| ), we get K\ = 2/rr\ 2 + 0(l/r 2 ), which, 
when substituted into the surface integral in (|6.6|), yields 



1 



1 



— / d^dtA-)dfAn ri = — + 0[- 
2n J\ x \<iz 13 \r 2 ) 13 r{ 2 \TZ 



As we can see, the integral is finite in the limit 1Z — > 00, with value 



(6.7) 



1 



The same method applied to Hi leads, since Hi = 0(l/r 2 ), to 



(6.8) 



2tt 



(6.9) 



We must compute now the spatial integral of 1/rf [see the first term in (|6.2| )1. It is clearly 
infinite because of the divergency at the bound x — > yi. By integrating 1/rf from r\ — e 
up to infinity, we obtain: J ri>£ <i 3 x/rf = 27r/e 2 , which is a pure constant fl57| , cancelled after 
applying the time-derivative in front of the 1/c-term in (|6.1|) . The second term in is 
therefore 

1 d ' ^ 3 ™ 2 ^ 



rf a x Iw£ VdV %V\ = -^^(ni 2 vi 2 ) + 0(2) + 1 ~ 2 . 

zcr 12 



(6.10) 



By summing the various contributions, we find the non-compact potential at 0.5PN order: 



£(NC) = _ G 3 m 2 m2 



1 



12r 3 



+ 



G mim 2 

cr\ 2 
G 3 m 2 m2 



2 J 

^12^12 



8r 2 r 2 2 ' 16 
1 
2 



+ ^-K 1 + H 1 



~lM^-5V)(n 12 vi 2 ) 



2cr 3 



ni 2 vi2) + 0(2) + 1 <->2 . 



(6.11) 
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Adding the other contributions in X we end up with the complete expression reported in 
the Appendix B. 

Finally, we give the value of the non-compact potential Q6.ll ) at the location of body 1. 



From the Hadamard recipe (|3.6|) we find for the functions K^ 2 and Hi >2 at point 1: 

2 



(Ki)i 
(#1)1 



3r? 



12 



3r? ' 



(#2)1 







12 



rf 2 ' 



(6.12a) 
(6.12b) 
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so 



that the non-compact potential at point 1 

G 3 m 2 



is 



(X( NC ))! = 



12 



1 2 1 2 

--m x + mim 2 + y^ m 2 



G 3 m\m 2 „, N 

H ^—3 — (-mi + m 2 )(ni2Ui2) + 0(2) 

zcr 12 



(6.13) 



VII. THE 2.5PN METRIC OF BINARY SYSTEMS 

From the results of the previous sections, we are in the position to write down the 
2.5PN harmonic-coordinate metric generated by two point-like particles as a function of 
the coordinate position x and a functional of the coordinate positions and velocities of the 
particles yi )2 (i), Vi )2 (t) (where t =const is the harmonic-coordinate slicing): 



0a«/CM) = ^AB/[x;yi(t),y 2 (*);vi(t), v 2 (t)] 



(7.1) 



The metric is given by (|2.5|) in which we insert the expressions of the potentials as listed in 
Appendix B. After combining together identical terms we obtain 1)8): 



fl'oo 



2Gm 1 
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+ G 2 m\m 2 



(-(n lVl ) 2 + Avf) 
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1 
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r 12 r 12 
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r 12 r 12 
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+ G ™T 2 f-|(ma«ia)^ - Q(n 12 v 12 )n\ 2 n{ 2 + 8ng,«2) + 0(6) + 1 <- 2 , (7.2c) 

(where we recall S = n + r 2 + ri 2 ). The harmonic-coordinate conditions satisfied by ([7.2[) 
read, with this order of approximation: ^[(— .g) 1 / 2 ^ ^] = 0(7) and <9 i/ [(— 5 , ) 1/ ' 2 fi ,li/ ] = 0(6). 

When applied to post-Newtonian initial conditions for the numerical evolution of two 
compact objects (provided that the initial spatial numerical grid does not extend outside 
the binary's near zone |4(|), the lengthy expressions ( |7.2j ) become a little simpler as the orbit 
can be considered as circular with a good approximation. In this case we have (^12^1) = 
0(5) = (rii2f2), the remainder 0(5) corresponding to radiation- reaction effects. Obviously, 
all the resulting 0(5)'s can be neglected since they yield terms falling into the uncontrolled 
remainders of ( |7.2| ). If in addition we are working in a mass centered frame, then y x = 
X 2 y\ 2 + 0(4) and y 2 = -Xiy 12 + 0(4), where X 1 = mi/m, X 2 = m 2 /m 1 m = m 1 + m 2 . 
All the remainders 0(4) become negligible after insertion in ( |7.2| ) . Thus, for instance: 
vj = Xlv 2 12 + 0(A), and 

(n 1 v 1 )=X 2 -(nv 12 ) + 0(A) , (7.3a) 

n 

(ni« 2 ) = -Xi-(rw;i 2 ) + 0(4) (7.3b) 

n 

(plus the same formulas with 1 and 2 exchanged). We denote by n = x/r and r = |x| 
the direction and the distance from the center of mass; r depends on the two individual 
distances r\, 2 through the relation 

r 2 = X ir \ + X 2 r\ - X x X 2 r\ 2 + 0(4) . (7.4) 

The magnitude of the relative velocity is v\ 2 = r i2 t^2PN + 0(6), where c^2PN denotes the 
orbital frequency of the circular motion at 2PN order, and is given by (8.6) below ||59|| . 

We provide also the values of the metric coefficients ( |7. 2|) computed at body 1 (since 
these might also be needed in the problem of binary coalescence), i.e. 

{gMt) = ^[yi(t);yi(t),y 2 (t);vi(t),v 2 (t)] , (7.5) 
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where the limit x — > yi is understood in the sense of ( |3.6|) . Directly from (7.2), or using the 
expressions of the potentials at 1 as given in Appendix B, we get: 
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Drastic simplifications occur in the case where the orbit is circular. 
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VIII. THE 2.5PN EQUATIONS OF MOTION OF BINARY SYSTEMS 



The motion of body 1 under the gravitational influence of body 2 is simply the geodesic 
motion taking place in the post-Newtonian space-time (|7.2j ). Now we write the geodesic 
equation of body 1 in the Newtonian-like form 



dt 



•1) 



The "linear momentum" vector and "gravitational force" (per unit of mass) are defined by 
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5.2) 



: 2 J 1 



The above quantities are computed using the regularization ( [Of ) which is a crucial ingredient 
of our point-mass model. 

Inserting the 2.5PN metric (|2.5| ) into (|8.2| ), we obtain 
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Next we replace into these expressions all the potentials and their gradients computed at 
point 1 as given in Appendix B (to this order the Hadamard partie finie is "distributive" 
), and get both V\ and T\ in terms of the relative separation y\ 2 = T\2n\ 2 and individual 



velocities v\ 2 [alternatively we can obtain V\ and T\ directly from ( |7.2|) 1. Then, we compute 
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the time-derivative of V[, and order-reduce all the resulting accelerations (which appear at 
orders 1PN or 2PN) by means of the 1.5PN equations of motion given by ( |3.11|) . After 
insertion in (|8 . 1|) and simplification, we end with the 2.5PN acceleration of body 1: 
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We find perfect agreement with the D amour- Deruelle |§-[F| equations of motion. To em- 
phasize the strength of this agreement we recall that the method employed in the present 
paper differs in many respects from the one originally used in @-|7] (see the discussion in 
the introduction). In the case of circular orbits, the equations reduce to 

dv\ 2 , , 32G 3 m 3 u 



~^2PN2/l2 



5c 5 rf- 



-v\ 2 + 0(6) . 
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dt UU"7 12 

The second term represents the standard damping force, while the orbital frequency uj 2 pn is 
the frequency of the exact circular motion at 2PN order, related to the harmonic-coordinate 
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separation r i2 by 



2 Gm 

a; 2PN = — 3~ 
'12 



1 + (-3 + i/) 7 + ^6 + y i/ + z/ 2 ) 7 2 



(8.6) 



Our notation is: m = mi + ms, v = mimj/m 2 (= X1X2 in the notation of the previous 
section) and 7 = Gm/rwc 2 . 



APPENDIX A: DERIVATION OF THE 2.5PN FLUID METRIC 



The derivation follows almost immediately from the results established in the section 



II. A of | 39fl . The Einstein field equations in harmonic coordinates are written as 



d v hT = , 
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where □ = ifd^dv is the flat d'Alembertian operator \rf v = diag(— 1, 1, 1, 1)], h^ u 
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fJ,U 



T]^ u , and A^ denotes the gravitational source term which is at least quadratic 



in h and its space-time derivatives (see for the expressions of the quadratic and cubic 
parts of A^ u ). From [^.^j, we have, to order 1PN, 
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Substituting the 1PN metric into the right-hand-side of the field equation we get 

\g\ = 1 + + 1 (w kk + 2V 2 ) + 0(6) , 

together with the gravitational source term (equations (2.12) in |3 

14 16 r 5 

A 00 = —jdiVdiV + ^{-Vd 2 t V - 2V i d t d i V + -(d t Vf 

+ \diVjidiVj + SdjVi) + diVdtVi - ^VdiVdiV 
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These are the needed equations, which lead, by application of the retarded integral on the 
right-hand-side of the field equations, to 
h oo + h u 2 „ 4 
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(the potentials are defined in the text). From this, we deduce the components of the covariant 
metric and find the result ( |2.5| ). It has been shown in section III of |39|] that the post- 
Newtonian metric matches in the external near zone to a solution extending up to the 
radiative zone. 



APPENDIX B: COMPLETE RESULTS FOR THE POTENTIALS 

We give first all the relevant potentials (valid all-over space-time) which are used in the 
obtention of the 2.5PN metric (ITT! 
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The values of the potentials at the location of body 1, following from (|3.6|) , are 
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The gradients of the potentials computed at body 1 (needed for the equations of motion) 
are 
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